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Abstract 

We consider the coherent exciton transport, modeled by continuous-time quan- 
tum walks, on Erdos-Reny graphs in the presence of a random distribution of 
traps. The role of trap concentration and of the substrate dilution is deepened 
showing that, at long times and for intermediate degree of dilution, the survival 
probability typically decays exponentially with a (average) decay rate which de- 
pends non monotonically on the graph connectivity; when the degree of dilution 
is either very low or very high, stationary states, not affected by traps, get more 
likely giving rise to a survival probability decaying to a finite value. Both these 
features constitute a qualitative difference with respect to the behavior found 
for classical walks. 
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1. Introduction 

Quantum walks provide a quantum extension of the ubiquitous classical ran- 
dom walks and have important applications in a broad range of fields including 
solid-state physics .polymer chemistry, biology, astronomy, mathematics and 
computer science [J 0, HI, 0] • Due to their coherent nature, the behavior of 
quantum walks can differ significantly from that of the classical random walks, 
as corroborated by measures of mixing times, hitting times and exit probabilities 
of quantum walks 

The continuous-time version of quantum walks (continuous-time quantum 
walks, CTQWs) has been extensively studied as effective model of energy trans- 
port in molecular systems such as chromophoric light-harvesting complexes 
[fil, 0|- I R photosynthesis, excitation energy is absorbed by pigments present 
in the antennas and subsequently transferred to a reaction center where an 
electron-transfer event initiates the process of biochemical energy conversion. 
This process has been studied for decades due to its impressive efficiency (even 
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over 90% in certain bacterial systems and higher plants), nonetheless a full de- 
scription of the mechanism leading to such a remarkable efficiency has not been 
achieved yet [8|, |9|. Also, the overall effect of the environment, of its (quenched) 
disorder and of the relative position of the reaction center are expected to play 
a central role ((j, Hoi ITU . Indeed, understanding how such a process works might 
be useful for the nano- engineering of optimized solar cells [9j . 

The success rate of an energy transfer process can be investigated by study- 
ing the interaction between a quantum walk (mimicking the rather coherent 
propagation of the exciton) with a reaction center (being it an impurity atom 
or molecule), which irreversibly traps the moving particle. Consequently, a great 
deal of recent theoretical work has focused on investigating essential features of 
basic trapping models, wherein a quantum particle moves in a medium contain- 
ing different arrangements of traps. The trapping problem on a one-dimensional 
structure has already been investigated in [12l . 1 1 31 ) , where distinct configurations 
of traps (ranging from periodical to random) where shown to yield strongly dif- 
ferent behaviors for the quantal mean survival probability, while classically, at 
long times, the exponential decay is always recovered. In this context the case of 
substrates displaying random topological inhomogeneity has not yet been inves- 



tigated, notwithstanding their experimental importance 14l. Il5j]. Such random 
structures, typically modeled by Erdos-Reny (ER) graphs, have attracted a 
great deal of interest in the last years also due to new tools introduced for their 
investigation (see e.g. [T6l 

In this work we study the survival probability of a CTQW moving on an 
ER graph in the presence of a fixed concentration of traps randomly placed. 
The number of traps M as well as the substrate connectivity encoded by the 
average number of neighbors per node z, are properly tuned in order to account 
for their role in affecting the trapping performance. Indeed, since ER networks 
lack hubs and display an overall homogeneous topology, the transport properties 
are controlled mainly by the average degree. Hence, for a given realization of 
the system we measure the survival probability Hm,z (t) as a function of time 
t, which, for sufficiently large systems, turns out to display a (qualitatively) 
robust behavior with respect to the realization of the system; however, due to 
the intrinsic randomness of the system (involving both the substrate topology 
and the trap arrangement) in order to properly outline the typical behavior, we 
generated several realizations (for given M, z and graph size N), over which we 
averaged to get (HM,z(t)}- Analogous calculations have been performed for the 
case of a classical particle modeled by a continuous-time random walk (CTRW) 
in order to figure out possible genuine quantum-mechanical effects. Analysis are 
performed both numerically and analytically relying on matrix diagonalization 
algorithms and on perturbation theory, respectively. 

Our results highlight that at long times and when the trap concentration is 
small, for large size (N » I) and intermediate degrees of dilution (z ~ O(N)), 
both quantal and classical survival probabilities typically decay exponentially 
with time. However, while the classical decay rate increases monotonically with 
z, the quantal decay rate displays a subtler dependence. Indeed, when the 
degree of dilution is either very large or very low, stationary states, not affected 
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by traps, get more likely and this makes the quantal survival probability to 
decay to a finite value related to the number of localized eigenmodes. As a 
result, given a fixed number of randomly arranged traps, in order to enhance the 
trapping efficiency we need to thicken the substrate connectivity, independently 
of the current degree of dilution, while quantum mechanically the strategy does 
depend on the current degree of dilution. 

This paper is organized as follows. In Sec. [2] we define CTRW and CTQW 
and we describe the substrate where they move on; then, in Sec. [3] we reports 
our results and finally in Sec. @] we discuss them and possible extensions. 

2. Coherent dynamics on random graphs 

The incoherent transport occurring over a discretizable environment can be 
modeled by continuous-time random walks (CTRWs) described mathematically 
by a master equation. However, when dealing with quantum particles at low 
densities and low temperatures, decoherence can be suppressed to a large ex- 
tent. Therefore, the study of transport in this regime requires abandoning the 
classical, master-equation-type formalism and adopt a quantum-mechanical ori- 
ented picture, where the local description of the complex network of molecules 
involved in the transport can be retained through a tight-binding approach. 

Interestingly, the CTRW picture can be mathematically reformulated to 
yield a quantum-mechanical Hamiltonian of tight-binding type; the procedure 
uses the mathematical analogies between time-evolution operators in statisti- 
cal and in quantum mechanics: The result are continuous-time quantum walks 
(CTQWs). 

In the following we provide the formal tools for the study of both CTRWs 
and CTQWs. 

2.1. Graph formalism 

Let us consider a graph Q made up of N nodes and algebraically described 
by the so-called adjacency matrix A: The non-diagonal elements A^ equal 1 if 
nodes i and j are connected by a bond and otherwise; the diagonal elements An 
are 0. We define the coordination number, or degree, of a node i as Zi — ^ ■ Ay. 

An Erdos-Renyi random graph is built as follows: Starting with N discon- 
nected nodes, every pair, say i and j, is connected, namely Aij = 1, with 
probability p, being < p < 1; multiple connections are forbidden and the 
extreme cases trivially correspond to a completely disconnected graph (p = 0) 
and to a fully connected graph (p — 1). In the limit of large size N, the coordi- 
nation number of an arbitrary node follows a binomial distribution with average 
z = p(N — 1) w pN. Due to their rigorous mathematical definition, ER graphs 
have been studied in details, from a topological point of view [HI, as well as for 
what concerns the properties of statistical mechanics models defined on them 
(see e.g. [3, HH). In particular, the Molloy-Reed criterion for percolation [2(| 
shows that, in the limit N — > oo a giant component exists, namely the graph is 
overpercolated, if and only if p is larger than 1/N. 
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Now, the Laplacian operator L of an arbitrary graph Q is defined as Ly = 
ZiSij — Aij, its spectrum, i.e. the set of all N eigenvalues of L, being denoted 
as = Ai < A2... < Aat; it follows from Gersgorin's theorem [2l|,|22| that L is 



positive semi-definite and, because its rows sum to 0, e^L = 0, where ejyr, is 
the row n-tuple each of whose entries is 1, therefore, the minimum eigenvalue 
is Ai = and it is afforded by eN- 

The eigenvalues and eigenvectors of Laplacian matrix L (also referred to 
as the admittance matrix, the stiffness matrix, or the Kirchhoff matrix) basi- 
cally form the backbone of any discussion of dynamic behavior of the networks 



represented by our graphs (see for example 23] and references therein). In 
particular, the degeneracy of the null eigenvalue represents the number of dis- 
connected components making up the whole graph (in the following we will 
always consider connected graphs made up of one single component), while the 
second smallest eigenvalue, also called spectral gap, controls the synchronization 
time. Moreover, the characterizations of eigenvectors are needed for a range of 
decentralized controls and dynamical- network analysis/design applications, in- 
cluding e.g., network partitioning, synchronization design, and optimal network 
resource allocation. Despite such need, graph theoretic studies of the Laplacian 



eigenvectors are sparse (see [24j for some reviews of these literature), and do 
not provide exact general characterizations of eigenvector-component values in 
terms of graph constructs for arbitrary graphs. 

In the following, we will refer to the eigenvalues and eigenvectors of L as the 
eigenvalues and eigenvectors of Q . 

To reveal the decay behavior of survival probabilities of the CTQW and 
CTRW we focus on systems of large size for which the spectral density of the 
Laplacian spectrum follows Wigner's law [25] and converges to the semicircle 
distribution 



where a — y Np(l — p) and z — pN. Moreover, for highly diluted (p <C 1) 
networks we can write a 2 ss z and p can be expressed as a function of z only. 

2.2. Classical and Quantum walks on graphs 

Continuous-time random walks (CTRWs) (2(| are described by the following 
Master Equation: 

d N 

being Pk,j{t) the conditional probability that the walker at time t is on node k 
when it started from node j at time 0. If the walk is unbiased the transmis- 
sion rates 7 are bond-independent and the transfer matrix T is related to the 
Laplacian operator through T = — 7L (in the following we set 7=1). 

We now define the quantum-mechanical analog of the CTRW, i.e. the 
CTQW, by identifying the Hamiltonian of the system with the classical transfer 
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matrix, H = — T 27J. Hence, given the orthonormal basis set \j), representing 
the walker localized at the node j, we can write 



N 



N 



H 



(3) 



j=l k~j 



where in the second term we sum over all connected couples k ~ j. The operator 
H is just the tight-binding Hamiltonian which applies to a large class of quantum 
transport systems such as excitons and charges in molecular and quantum dots 



Now, the quantum mechanical time evolution operator is defined as U(£, to) = 
exp[— iH(t — to)], so that the transition amplitude a k ,j{t) from state \ j) at time 
to state |fc) at time t reads ct k ,j(t) = (k\U(t,0)\j), and it obeys the following 
Schrodinger equation: 



dt 



a k,j(t) 



N 

' E "if '/ ; :/ : - 

1=1 



(4) 



formally very similar to Eq. [2] Then, the classical and quantum transition 
probabilities to go from state \ j) to state in a time t are given by Pk j(t) = 
(fc|e- tT |j) and n k<j (t) = \a kj (t)\ 2 = \(k\ e - itH \j)\ 2 , respectively. 

In the absence of traps and other impurities, the operators describing the 
dynamics of CTQWs and of CTRWs share the same set of eigenvalues and 
of eigenstates; denoting with E n and |$ n ),n € [1,N] the nth eigenvalue and 
orthonormal eigenvector of L, we can write 



At)=J2e- E « t {k\3> n ){$ n \ 



(5) 



and 



n k,j(t) 



N 



(6) 



2.3. CTRWs and CTQWs in the presence of traps 

Let us introduce a set M of M traps placed randomly on nodes {mi, TO2, wm}- 
In the incoherent, classical transport case trapping is incorporated into the 
CTRW according to 

T = T -r = -L-r, (7) 

where To denotes the unperturbed operator without traps while T is the trap- 
ping operator defined as 



\m) (m\ 



(8) 



The capture strength T determines the rate of decay for a particle located at 
trap site and here it is assumed to be equal for all traps. 
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The transfer operator T is therefore self-adjoint and negative definite; we 
denote its eigenvalues by —A; and the corresponding eigenstates by \4>i). 
The mean survival probability for the CTRW can be written as 



P M (t) = 



1 



E 



N-M ^ ^ 



Pkj(t) 



1 



N-M 



N 

E< 

1=1 



E w< 



k4M 



(9) 



From Eq. [9] one may deduce that Pm (t) attains in general rather quickly an 
exponential form; furthermore, if the smallest eigenvalue A m i n is well separated 
from the next closest eigenvalue, PmU) is dominated by A m i n and by the corre- 
sponding eigenstate |0 m i n ) fl2lfl3j|: 



P M (t) 



1 



N-M 



-A min t 



E w« 



(10) 



As for quantum transport, in the presence of substitutional traps the system 
can be described by the following effective (but non-Hermitian) Hamiltonian 



12] 



H = Hi 



iT, 



(11) 



where Ho denotes the unperturbed operator without traps. 

Due to the non-hermiticity of H, its eigenvalues are complex and can be 
written as E[ — e/ — iji (I — 1,...,N); moreover, the set of its left and right 
eigenvectors, |$/) and ($;|, respectively, can be chosen to be biorthonormal 



((<£>z|<I>J) = Siji) and to satisfy the completeness relation Ylt=i = !• 

Therefore, according to Eq. 01 the transition amplitude can be evaluated as 



N 

E' 

i=i 



-ill +ie ( )i 



(12) 



from which iVk,j[t) — \&k,j(t)\ 2 follows. 

Of particular interest, due to its relation to experimental observables, is the 
mean survival probability n^i) which can be expressed as [l2| 



n 



M 



N ■ 



m E E 

j$M k£M 
N 



7Tfcj(t) 



N-M ^ 

i=i 



N 

M E ' 

l,V = l 



2 2^ ($,|m)<m|*j) 



-i{Ei-E"[,) 



N 




$//|m)(m|$ z ) 



(13) 
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The temporal decay of IhAt) is determined by the imaginary parts of Ei, i.e. 
by the 7;. As shown in [12], at intermediate and long times and for M <C N, 
the Tlm (t) can be approximated by a sum of exponentially decaying terms: 



1 N 



.^E'-'. (") 

and is dominated asymptotically by the smallest 7; values. 



3. Results 

As shown in the previous section, the long-time decay of the survival prob- 
abilities PmW an d n M (i) is controlled by the "spectra" {Aj} and {7/}, respec- 
tively. When the capture strength is small (r -C 1), some insights into such 
spectra can be obtained following a perturbative approach to get the correction 
to the unperturbed eigenvalues. For instance, the first-order correction to Ex, 
corresponding to the (unperturbed) eigenvector |4>j °') = e^/v^/V, reads as 

V- l/„i*(0)\ 2 - M 



E^=-iT £ |H$i 0) ) = -* r 77F' (15) 



where, recalling that the graph is connected, we applied the perturbative theory 
for non-degenerate eigenvalues. 

Let us first focus on the classical case. From Eq. [15] one can also write 

A 1 = r^ + o(r 2 ), (16) 

and, similarly, for the remaining eig envalues A, = A, (0) + rA, (1) + 0{T 2 ). Now, 
for large system size the smallest non-null eigenvalue \% can be estimated as 
A^ ~ z — 2a — Np — 2yjNp{l — p) (see Eq. [T]), which for large average degree 
2, is well separated by Ai, even if corrected by a term order of T. Therefore, 
in this case the survival probability is controlled by A m ; n = Ai ~ TM/N. On 
the other hand, when the dilution gets lower the spectral gap also decreases 
and the smallest perturbed eigenvalue A m i n can get smaller than TM/N. Of 
course, the lower z and the smaller the expected A m j n , in such a way that 
the survival probability decreases with a smaller rate. This is confirmed by 
numerical calculations: As shown in Fig. [T] the numerical data pertaining to a 
given realization of the substrate are well fitted by the function 

PmM ~ exp(-r© c t), (17) 

moreover the exponent O c gets lower when the link probability is drastically 
reduced (see Fig. [3]). 

The quantum-mechanical case is more subtle as in general the first eigenstate 
\^^) is not sufficient to determine the smallest perturbed eigenvalue of the 
complex spectrum. 



7 



Figure 1: Average Classical and Quantum survival probability as a function of time on a 
semilogarithmic scale-plot; while N = 40, M = 1 and T = 0.1 are kept fixed, different values 
of p are considered, that is, p = 0.11 (darker color for the quantal probability and «, for the 
classical probability), p = 0.5 (intermediate color for the quantal probability and □, for the 
classical probability), and p = 0.95 (brigheter color for the quantal probability and A, for 
the classical probability). The dashed, clear line represents Eq. 1171 with C = M/N and is in 
full agreement with numerical data. In the inset we show only 11m, 20 M on a semilogarithmic 
scale plot in order to highlight the exponential decay; the dashed line represents the best fit 
with O q ~ 10 -5 . These data refer to a particular realization of the substrate. 



First of all, we notice that graphs displaying either high or low link dilution 
are more likely to present (at least one) Laplacian eigenvector displaying a large 
number of null entries. This is easy to see by recalling, respectively the "Edge 
Principle" jU Q and the fact that for a complete graph K jv the eigenvectors 
show a monotonic trend in confinement, that is, as we go toward faster modes, 
the eigenvector is more localized and the region of support decreases by one 
node in the graph. 

Another, intuitive, way to see this point is by noticing that the so-called Faria 
vectors are more likely to be Laplacian eigenvectors in the above mentioned 
regimes. Indeed, we recall that, in the "valuation notation" [28j . a Faria vector 
|£) is a vector with nonzero entries only on two vertices i and j with = 
— = 1; moreover, a Faria vector is an eigenvector of the Laplacian of the 
graph Q if and only if i and j are twins, i.e., if every vertex v (fc {i,j} is either 
adjacent to both i and j or to neither one of them. The corresponding eigenvalue 
is A = Zi + 1 = Zj + 1 if G E(Q) and \ = Zi = Zj if [i, j) £ E(Q). Now, the 
probability that such conditions are fulfilled can be written as 




k=l 



(1-P) 



2N-4 



+ [1- 2(1 -p)p] 



N-2 



(18) 



which is not negligible only in the region of very high and very low dilution and 
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Figure 2: Comparison between the mean average classical and quantum survival probability 
(Pm.z) and (TIm respectively, for a system made up of V = 40 nodes in the presence of 
a random distribution of M = 4 (upper panel) and M = 16 (lower panel) traps. Several 
degrees of dilutions are considered as shown in the legend, common for both panels. Notice 
that curves corresponding to the classical case are overlapped. Averages have been performed 
over 10 3 different realizations. 



for relatively small sizes; in particular, when p scales like p = j/N, being 7 a 
finite value, from Eq. [T5]we get Pp ~ exp(— 2j)j 2 /N for N 1. Hence, a Faria 
vector is more likely to be an eigenvector as p approaches either 1 or 0, namely 
when there exists relative (local) homogeneity for the two nodes corresponding 
to the non-null entries. 

According to formula in Eq. ll5l the first-order correction is zero whenever the 
traps are positioned in any node other than the two twin sites; analogously, one 
can verify that higher order corrections are null, indeed, such highly localized 
states do not see the traps at all. Consequently, the average survival probability 
decays to a finite value, similarly to what evidenced in 13[ when deterministic 
and regular arrangements of traps were considered. The asymptotic value of 
(J^-m,e) is given by the normalized number of stationary modes not intercepting 
the traps which, neglecting higher order terms, can be estimated by the number 
of disjoint Faria twins uf times the probability that they are not occupied by 
traps, namely (1 - M/N) 2nF 0. 

The above picture is confirmed by numerical results shown in Fig. [TJ where 



1 Nonetheless, a rigorous estimate should take into account that high dilution can induce 
the disconnection of a subset of nodes, while at low dilution correlations among twins cannot 
be neglected. 
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Figure 3: Ratio between the mean average classical and quantum survival probabilities (Pm,z) 
and (JIm z)j pertaining to different degree of dilutions, z and z' , for a system with V = 400 
and M = 1. We choose as reference value z = N/2, while z' is taken either very small or 
comparable with the system size, as shown in the legend. Notice that the case p = 0.5 gives 
rise to relatively small survival probability in the quantal case, while classically it corresponds 
to an intermediate trapping efficiency in agreement with intuition. 



the qualitative difference with respect to the classical case can also be noticed. 
On the other hand, for large sizes and intermediate degree of dilution, all the 
simulations performed recover the expected exponential decay at long times (see 
Eq . [Til) . Indeed, numerical data can be properly fitted by the function 

n A f(<)«exp(~re 9 t). (19) 

Finally, we performed averages over several realizations of the underlying 
graph, so to get the mean average probabilities (Pm,z) and (Hm,z)- 

Results are shown in Fig. [5] and Fig. [3J in the latter figure we also plotted 
the ratios (P M ,z{t)) / '<*Vs'(*)> and ^mA*)) / ^M.g'ff)). B Y fi xing z = N/2, 
hence corresponding to an intermediate dilution, and z' <C N or z' ~ N, we see 
that, classically, by reducing the dilution, the survival probability gets larger 
and vice versa, as expected. On the other hand, in the quantal case, both high 
and low dilution regimes prove to be less efficient in trapping the particles. 
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The slow decay we highlighted for low- and high-connectivity networks is 
consistent with results found previously [29| for the long time average x of a 
CTQWs embedded in ER random graphs. It is worth recalling that x is defined 
as 



where is the long time averaged transition probability which, classically, 
equals to the equal-partionend probability 1/N, while quantum-mechanically is 
given by 



Now, x is found to be larger than 1/N and to be almost a constant value in a 
wide range of average degree z but increases slightly as p approaches from above 
the percolation threshold and it increases fast when the network approaches 
a fully connected one where x = 1 — 0(N~ 1 ) j29j. Otherwise stated, high 
(and smaller) connectivities correspond to a large degree of localization hence 
reducing the probability to get trapped. 

4. Conclusions 

In this work we considered a continuous-time quantum walk (CTQW) prop- 
agating on Erdos-Renyi random graphs of size N endowed with a tunable link 
probability p, in such a way that the average coordination number is z = pN; 
moreover, M sites extracted randomly are occupied by traps. We measured 
the survival probability HM,z{t) and we compared it to the analogous survival 
probability Pm,z (t) found for a classical continuous-time random walk (CTRW). 
As expected from analytical arguments, when M <C N and the graph displays 
a relatively large size with intermediate degree of dilution, in the long-time 
regime both functions typically decay exponentially with time. However, when 
the dilution degree is either very low or very large, while Pm,z(£) still decays 
exponentially, IIif ]Z -(t) can decay to a finite value due to the existence of local- 
ized eigenmodes not "perceiving" the traps; of course, this effect gets less likely 
as M is increased. 

We stress that the crossover evidenced in the behavior of Hm,z for different 
regions of z stems from a different degree of localization exhibited by the quan- 
tum walk. Indeed, in networks with either large or very small mean degree, 
the quantum excitement is on average most likely to be found at the initial 
node |3(3] • The reason is that the substrate topology allows the establish- 
ment of eigenstates supported by a restricted subset of nodes; this is of course 
a quantum- mechanical effect: the classical particle quickly reaches an equipar- 
tition condition being equally spread on the whole substrate. 

As a result, while the efficiency (in terms of likelihood of trapping) of classical 
transport is rather robust with respect to the degree of dilution but it can be 
substantially reduced by cutting links, for quantum transport, improving the 




(20) 




(21) 
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efficiency by adding/removing links is more a subtle operation, as it sensitively 
depends on the starting topology. 
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